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Abstract 

Bees and wasps could exhibit shape and size sexual dimorphism, and most of their 
morphological variation could depend on phenotypic responses due to environmental pressure 
during ontogenetic development. More complex measurement techniques related to size and 
shape rather than simply to mass and length should be required to analyze such a complex sexual 
dimorphism. In this study, differences related to wing shape and size of males and females of 
Ophion intricatus Brulle (Hymenoptera: Ichneumonidae) were evaluated using geometric 
morphometries. Right and left wings of specimens were used, and a photographic matrix was 
constructed in which 18 morphological landmarks based on shape and vein patterns of the wings 
were digitalized. A multivariate analysis of wing shape showed significant differences between 
sexes and sites. The geometric variation demonstrated that the points at the intersection of radial 
and cubital-anal veins might be key characters to differentiate between sexes. This study also 
showed the presence of two clearly different male morphotypes coexisting in the same study site. 
However, it should be noted that the results of this study showed that the variation in wing shape 
is an analytical character in the determination of sexual differences in the family Ichneumonidae. 
These differences raise the question of whether sexual dimorphism of wing shape may be 
modulated by natural selection. 
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Introduction 



Bees and wasps, as most insects, exhibit 
sexual size and shape dimorphism (e.g., 
Hurlbutt 1987; Mackauer 1996; O'Neill 
2001; Pretorius 2005; Mitrovski-Bogdanovic 
et al. 2009; Zikic et al. 2009). Most of the 
morphological variation in wasps could be 
due to phenotypic responses (plasticity) that 
act during ontogenetic development as a 
consequence of environmental pressures 
(Shreeves and Field 2008). Females have 
some adaptive advantages over males, such 
as greater fecundity and parental care, be- 
cause females are typically larger than males 
(Forrest 1987; Andersson 1994; Moller and 
Zamora-Munoz 1997; Reeve and Farbain 
1999). However, in some species, males are 
longer but have less mass than females (e.g., 
Cepeda-Pizarro et al. 1996). Those consider- 
ations suggest that inferring and quantifying 
sexual dimorphism requires more complex 
measurement techniques related to shape and 
size rather than simply to mass and length 
(Gidaszewski et al. 2009; Benitez et al. 
2010). Another important factor to be evalu- 
ated within sexual dimorphism is allometry, 
defined as the association between size and 
shape, or the covariation of parts due to vari- 
ation in size. All types of allometry can 
contribute to the sexual differentiation at the 
respective levels of variation (Gidaszewski 
et al. 2009). Allometry is classified accord- 
ing to the cause of variation in size that gives 



rise to the allometric relationship, which is 
usually ontogenetic growth or evolutionary 
change; therefore, there is either ontogenetic 
or evolutionary allometry, respectively. If 
individuals of the same species are examined 
at the same ontogenetic stage, both of these 
sources of "dynamic" variation are kept con- 
stant, and the remaining individual variation 
would be "static," corresponding to static 
allometry (Klingenberg 1996). 

Ophionines are solitary koinobiont endopar- 
asitoids of the larvae of other 
holometabolous insects. Many ophionines 
oviposit into relatively mature larvae in their 
3 rd or 4 th (Vickery 1929) or 4 th or 5 th instar 
(Moutia and Courtois 1952; Yu and Horst- 
mann 1997), though some species attack 
young larvae (Price 1975). Generally a sin- 
gle egg is deposited free in the haemocoel. 
No species of ophionine are known to devel- 
op gregariously (Gauld and Bolton 1988; Yu 
and Horstmann 1997). The genus Ophion 
Fabricius (Hymenoptera: Ichneumonidae) 
includes about 200 species (Gauld 1978) dis- 
tributed worldwide, and most of them 
inhabit temperate climates (Townes 1971). 
Most Ophion species are parasitoids, laying 
their eggs in Lasiocampidae, Noctuidae, and 
Notodontidae pupae or larvae. 

In the present study, static allometry is ana- 
lysed on the species Ophion intricatus Brulle 
(Hymenoptera: Ichneumonidae). Geometric 
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morphometries techniques were used, with 
the aim to investigate sexual dimorphism of 
wings and to characterize the wings as a trait 
of sexual differentiation, as well as to show 
that there is a differentiation pattern in males 
of the same species coexisting in the same 
habitat. This is the first step for further stud- 
ies on behavior at group level in the area. 

Materials and Methods 

Data acquisition 

A total of 136 specimens of O. intricatus (50 
female and 86 males) were used in this 
study, and they were collected in the Botani- 
cal Park "Pedro del Rio Zanartu" - Bio Bio 
Region, Chile (36° 47' S, 73° 09' W). Wasps 
were collected with UV light traps (GE 
Power Pro 950, www.ge.com ) using an 800- 
watt electric generator. The light sources 
were placed over a white sheet to increase 
the luminosity. Traps were installed for 4 hr 
in different sampling points. Following the 
key of Gauld (1978), specimens were care- 
fully checked for their species status, and the 
process was based on taxo-morphological 
characters. Once collected, specimens were 
processed and mounted. Right and left wings 
were removed from the body and then 
mounted in fixed molds (compressed be- 
tween two microscope slides) for subsequent 
image capturing. 
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(Zelditch et al. 2012). For this purpose, 
TpsDig 2.12 software (Rohlf 2008b) was 
used. The x-y coordinates of the biologically 
homologous landmarks were then aligned 
and superimposed with the method of mini- 
mum least squares based on the generalized 
procrustes analysis (Dryden and Mardia 
1998; Rohlf and Slice 1990). 

Once the Cartesian x-y coordinates for all 
points were obtained, shape information was 
extracted from the coordinate data by using a 
full Procrustes fit (Rohlf and Slice 1990; 
Dryden and Mardia 1998). Due to the 
matching symmetry of the wings, the dataset 
contained two separate landmark configura- 
tions, one for the left and one for the right 
side. The analysis of shape included the re- 
flection of all configurations from one body 
side to its mirror image (for details see 
Klingenberg and Mclntyre 1998; Klingen- 
berg et al. 2002). After reflection, all 
configurations can be superimposed in the 
Procrustes fit, a procedure that removes in 
formation on location orientation and rota- 
tion and standardizes each specimen to unit 
centroid size, which is the square root of the 



Table I. Anatomical description of landmarks in the wing 
of Ophion intricatus. A = Anal; Cu = Cubital; M = Medium; R 
= Radial; Rs = Sector radial; M+Cu = Medium cubital; m = 
Medium transversal; cu = Cubital transversal; m+cu = 
Medium cubital transversal. 



Landmark 


Anatomical structure 


#1 


Proximal Base M+Cul 


#2 


Distal Base 1A 


#3 


Distal Base M+Cul 


#4 


Distal Base A 


#5 


Proximal Base Culb 


#6 


Proximal Base lm-cu 


#7 


Prior Base M 


#8 


Prior Base of stigma 


#9 


Distal Base Culb 


#10 


Proximal Base M 


#11 


Prior Base Rs 


#12 


Proximal Base M 


#13 


Distal Base Rs 


#14 


Distal Base Cu2b 


#15 


Distal Base Cula 


#16 


Distal Base M 


#17 


Prior Base Rs 


#18 


Proximal Base Rl 



Morphometric analysis 

Right and left wings of wasps were photo- 
graphed with a Sony 10 DSC-H7 camera 
( www.sony.com ) with directed fiber optics 
lighting. Then, photographic matrixes were 
constructed using TpsUtil 1.57 program 
(Rohlf 2008). Eighteen morphological land- 
marks (Type 1 landmarks) were digitized for 
left and right wings in all specimens, select- 
ed according to external anatomy and vein 
pattern of the wings (Figure 1, Table 1) 
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summed squared Euclidean distances from 
each landmark to the specimen centroid and 
provides an estimate of the size of the stud- 
ied structure (Dry den and Mardia 1998). 
Centroid size was used as a measure of 
whole body size, while the information 
about the shape variation was extracted from 
the Procrustes superimposition. 

After performing the Procrustes fit, the co- 
variance matrix generated from the 
Procrustes coordinates, and which included 
the measures of the association between 
Procrustes coordinates (the x-y coordinates 
of each landmark after the Procrustes super- 
imposition), were extracted. Shape variation 
in the entire dataset was then analyzed with 
principal component analysis (PCA) based 
on decomposition of covariance matrix. 

Additionally, to assess shape and size varia- 
tion between sexes, Procrustes ANOVA was 
performed, and the results are reported as 
sums of squares and mean squares that are 
dimensionless (Klingenberg and Mclntyre 
1998; Klingenberg et al. 2002). 

Static allometry was evaluated by a multi- 
variate regression of shape, pooled within 
sex (Procrustes coordinates), on centroid 
size. Then, the predicted values were used 
as shape variables, accounting for the al- 
lometric component of shape variation, and 
the residuals were used as the non-allometric 
component of shape variation. A plot of the 
regression scores (the projection of shape 
data onto the direction of the regression vec- 
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tor in the shape tangent space) was per- 
formed to better visualize allometric rela- 
relationships (Drake and Klingenberg 2008; 
Klingenberg et al. 2012). All morphometric 
and statistical analyses were performed us- 
ing the free morphometric software Morpho 
J 1.05f (Klingenberg 2011). 

Results 

PCA for allometric and non-allometric com- 
ponents of within-sex variation, showed how 
the first principal components accumulated 
shape variation in very few dimensions 
(Figure 2). 

Regarding allometric components, the first 
three PCs accounted for 46.6% (PCI + PC2 
+ PC3 = 18.178% + 16.365% + 12.061%) of 
the total shape variation, and provided a rea- 
sonable estimation of the total variation. The 
other PCs each accounted for no more than 
8.8% of the variation. 

Regarding non-allometric components, PCI 
held more than 30% of the complete varia- 
tion between sexes, and the remnant PCs 
each accounted for no more than 13.6% 
(Figure 3, 4). Procrustes ANOVA showed 
that variation in size and symmetric shape 
between sexes was highly significant (Table 
2). To quantify static allometry, multivariate 
regression of shape (pooled within sexes) 
was used on centroid size. The results con- 
firmed the occurrence of a strong static 
allometric effect that accounted for 30.48% 
of the total shape variance (p < 0.0001). 



Table 2. Procrustes ANOVA for both centroid size (CS) and shape (SH) of Ophion intricatus. Sums of squares (SS) and mean 
squares (MS) are in units of Procrustes distances (dimensionless). 





SS 


MS 


df 


F 


P 


Pillai tr. 


p (param) 




Sex 


0.000051 


0.000025 


2 


133.5 


<0.0001 






CS 


Individual 


0.000012 


0.0000001 


65 


2.09 


0.0016 






Side 


0.000001 


0.000001 


1 


1.55 


0.217 








Ind*Side 


0.000006 


0.0000001 


67 












Sex 


0.11457622 


0.001790254 


64 


42.37 


<0.0001 


1.68 


O.0001 


SH 


Individual 


0.08788273 


4.22513E-05 


2080 


3.5 


<0.0001 


22.16 


O.0001 


Side 


0.00251535 


7.86046E-05 


32 


6.51 


<0.0001 


0.83 


O.0001 




lnd*Side 


0.02587538 


1.20687E-05 


2144 
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Discussion 

The results of this study confirmed the pres- 
ence of sexual variation in wing shape in O. 
intricatus. A number of studies have found 
that sexual dimorphism associated with size 
and shape in insects is adaptive (Nylin and 
Gotthard 1998; Farbain et al. 2007; 
Gidaszewski et al. 2009; Benitez et al. 
2013). Some of the hypotheses proposed to 
explain this kind of dimorphism in arthro- 
pods do not take into account whether the 
differences could be due exclusively to sex- 
ual dimorphism or to allometric processes 
(Walker and Rypstra 2001; Esperk et al. 
2007; Gidaszewski et al. 2009). According 
to Zikic et al. (2010), allometry contributes 
to the variation in body shape among female 
biotypes of Ephedrus persicae. Several mor- 
phological characters, such as size of head, 
width of petiolus, width of mesoscutum, and 
length of ovipositor sheath, diverged among 
the E. persicae biotypes. In our study, the 
influence of size on sexual differences was 
taken into account, and also if there was 
some effect on the shape differentiation 
(Klingenberg 1996; Fairbairn 1997; Debat et 
al. 2003). 

The wings of O. intricatus were evaluated, 
and it was determined that they were strong- 
ly under sexual dimorphism. The results 
show that females tended to have longer 
wings than males, as seen in the displace- 
ment of landmarks 4, 7, and 8 (the 
connection of the transversal radial and cubi- 
tal-anal veins). This type of wing 
dimorphism is present in almost all parasi- 
toids within Ichneumonidae, as the females' 
metasomas are likely to be a bit heavier than 
males because the females carry eggs. This 
extra weight could result in slight compensa- 
tory differences in wing size and shape. This 
kind of variation may be an adaptive re- 
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sponse to flight, and for behavioral activities 
such as localization, foraging, and oviposi- 
tion. In contrast, males' wings are wider, 
which might be linked to the selection of 
courting sites (longer flight periods in search 
of sexually active females) and to predator 
avoidance. 

A differentiation of wing shape between two 
male populations in the same sampling site 
was found, and two male morphotypes 
seemed to coexist. The presence of these two 
male morphotypes might be the effect of 
sexual differentiation due to females not pre- 
ferring a particular male wing shape. 
However, it must be taken into account that 
male populations used different microhabi- 
tats. The light traps used could reasonably 
have attracted specimens from different mi- 
crohabitats during the collecting process, 
and males may have morphologic micro- 
adaptation to their habitat. Zikic et al. (2010) 
reported that wing shape characters related 
to allometry differed among female speci- 
mens of E. persicae that emerged from 
different aphid hosts, indicating differences 
in their ecology and physiology. 

It is expected that the anatomical design of 
the wings of insects will be subjected to 
strong selection. Here, several wing mor- 
photypes emerged, raising the question of 
whether sexual dimorphism on wing shape 
may be modulated by natural selection 
(Betts and Wootton 1988; Dudley 2000). 

Behavioral studies are needed in order to 
clarify what the adaptive significance of the 
findings of our study are. However, it is 
worth noting that the results of this study 
showed that variation in wing shape was an 
analytical character in the determination of 
sexual differences in the individuals of O. 
intricatus. 
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Allometric 



T I r 



Principal components 



Non-allometric 



1 



10 15 20 25 



Principal components 

Figure 2. Total shape variation by PCA is shown, using covariance matrices of a) individual variation (allometric component) 
and b) corrected by size (non-allometric component). High quality figures are available online. 
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Figure 3. PCA of wing shape variation (allometric component, pooled within sexes). The different confident ellipses show fe- 
males and two male morphotypes with a 95% level of confidence. In the x and y axes, PCI and PC2 are respectively shown with 
associated shape deformation images. High quality figures are available online. 




Figure 4. PCA of wing shape variation (non-allometric component, pooled within sexes). The different confident ellipses show 
the female and two male morphotypes with a 95% level of confidence. The first two PC axes with shape deformation images 
associated are shown. High quality figures are available online. 

V ) 
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